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I. INTRODUCTION 



Considerable attention has been focussed on the infrared fixed point behavior 0| of the standard model (SM) and 
its extensions, especially the minimal supersymmetric standard model (MSSM). This is because there may be a stage 
of unification beyond the SM, and if so, it then becomes important to perform the radiative corrections in determining 
all the dimension < 4 terms in the lagrangian. This can be achieved by using the renormalization group equations 
to find the values of parameters at the electroweak scale, given their values at the unification scale. As such, much 
effort has been devoted to the study of the evolution of various dimensionless Yukawa couplings in the SM, and its 
minimal supersymmetric extension, the MSSM. Using the renormalization group evolution, one can relate the Yukawa 
couplings to the gauge couplings via the Pendleton- Ross infrared stable fixed point (IRSFP) for the top-quark Yukawa 
coupling [^js), or via the quasi-fixed-point behavior Q|. The predictive power of the SM and its supersymmetric 
extensions may, thus, be enhanced if the renormalization group (RG) running of parameters is dominated by infrared 
stable fixed points (IRSFPs). These parameters (Yukawa couplings, ratios of Yukawa couplings to gauge couplings, 
etc.) do not attain fixed point values at the weak scale, the range between the grand unified theory (GUT) scale 
and the weak scale being too small for them to closely approach the fixed point. Nevertheless, the couplings may be 
determined by the quasi-fixed point behavior Q where the value of the coupling at the weak scale is independent of 
its value at the GUT scale, provided the coupling at the GUT scale is large. 

In supersymmetric theories there are superpartners of the ordinary particles in the spectrum, due to which there 
are additional Yukawa couplings Q in these models which lead to baryon number (B) or lepton number (L) violation. 
Often, a discrete symmetry Q called i?-parity (Rp) is introduced to eliminate these B and L violating Yukawa 
couplings. However, the assumption of Rp conservation in MSSM appears to be ad hoc, since it is not required for 
the internal consistency of the model. Considerable attention has, thus, recently been focussed on the study of the 
renormalization gro up evolution of Yukawa couplings, including the baryon and lepton number violating couplings, 
of the MSSM 0^]. This includes the study of quasi-fixed-point behavior as well as the true infrared fixed points 
of the different Yukawa couplings, and the analysis of their stability. It has been shown that in the Yukawa sector 
of the minimal supersymmetric standard model with baryon and lepton number violation there is only one infrared 
stable fixed point. This corresponds to nontrivial fixed point for the top- and bottom-quark Yukawa couplings and 
the B violating coupling A233 , and a trivial one for the r- Yukawa coupling and the L violating coupling A233 . It was 
shown that all other fixed points are either unphysical or unstable in the infrared region. Similarly, fixed points were 
obtained for the corresponding soft supersymmetry breaking trilinear couplings as well [^. 

The purpose of the present paper is twofold. Firstly, we study the renormalization group evolution in the minimal 
supersymmetric standard model with baryon and lepton number violation in order to obtain the exact as well as 
the approximate analytical solutions for the Yukawa couplings and the soft supersymmetry breaking couplings at the 
weak scale given their initial values at the ultraviolet (UV) or the GUT scale. Second, we study the renormalization 
group flow of such a system, and determine the infrared fixed surfaces and infrared fixed points toward which the RG 
fiow is attracted. 

The plan of the paper is as follows. In Sec. II we describe the renormalization group equations for the minimal 
supersymmetric standard model with baryon and lepton mumber violation involving the highest generations. We 
obtain the exact solutions for the RG equations, and describe the infrared fixed points for the third generation 
Yukawa couplings and the highest generation baryon and lepton number violating couplings. Here we also study the 
corresponding RG equations for the soft supersymmetry breaking trilinear couplings, and describe their infrared fixed 
points. Sec. HI is devoted to the study of the approximate analytical solutions of the RG equations for the Yukawa 
couplings and the soft supersymmetry breaking trilinear couplings. In Sec. IV we carry out a detailed numerical study 
of the infrared attractive fixed surfaces, and present their two dimensional projections to demonstrate the existence 
of a strongly attractive fixed point. In Sec.V we present the summary and conclusions. 



II. RENORMALIZATION GROUP EQUATIONS AND INFRARED FIXED POINTS 

A. Infrared fixed points for Yulcawa couplings 

In this section we recall some of the basic features of the renormalization group evolution in the minimal supersym- 
metric standard model with baryon and lepton number violation, and obtain the infrared fixed points of the Yukawa 
couplings and the soft supersymmetry breaking trilinear couplings of the model. The superpotential of the model is 
written as 

W = {hu)abQlURH2 + {hD)abQl^RHi + {hE)abLlE\Hi + fiHiH2, (1) 
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where L, Q, E, D, U denote the lepton and quark doublets, and anti-lepon singlet, d-type anti-quark singlet and 
u-type anti-quark singlet, respectively. In Eq. (l]), {hu)ab, {hD)ab and {hE)ab are the Yukawa coupling matrices, with 
a, 6, c as the generation indices. Gauge invariance, supersymmetry and renormalizability allow the addition of the 
following L and B violating terms to the MSSM superpotential (0): 



Wl — -XabcLIL'lEjj + X'abcLIQ'iD + fiiLiH2, 
Wb = iKbc^RDnUR, 



(2) 
(3) 



respectively. In this paper we shall consider, apart from the Yukawa couplings in (|1|), the dimensionless Yukawa 
couplings Xabc, X'^bc ^^"^ '^'abc o^ily, and ignore the dimcnsionful couplings fj, and /i,;. The couplings Xabc and A^'^^ are 
antisymmetric in their first two indices due to SU{2)l and SU{3)c group structures, respectively. Corresponding to 
the terms in the superpotentials (|l|) , (Q) and , there are soft supersymmetry breaking trilinear terms which can be 
written as 



Vsoft = {Au)ab{hu)abQlUjiH2 + {AD)abihD)abQlDjjHi 



+ {A-E)abihE)abL1Ej^Hi 



■^{A\)abcXabcLlL\E If + {Ax')abcX'abcL'LQ)^D 

1 ^ b ~ c 

2 {A\" )abcKbcDRDRU R 



(4) 



where a tilde denotes the scalar component of the chiral superfield, and the notation for the scalar component of 
the Higgs superfield is the same as that of the corresponding superfield. In addition there are soft supersymmetry 
breaking gaugino mass terms with the masses Mi with i = 1, 2, 3, corresponding to the gauge groups SU{2)l, 
and SU{3)c, respectively. 

Since the third generation Yukawa couplings are the dominant couplings in the superpotential (|^) , we shall retain 
only the elements (/i£/)33 = /it, (/id)33 = hb, (/il)33 = h-r in each of the Yukawa couplings matrices hu, hu, h^, 
setting all other elements equal to zero. Furthermore, since there are 36 independent L violating trilinear couplings 
Xabc and AJjf,^ in (||), and 9 independent B violating couplings A"^^ in the baryon number violating superpotential 
(^, we would have to consider 39 coupled nonlinear evolution equations for the L violating case and 12 coupled 
nonlinear equations for the B violating case, respectively. Thus, there is a clear need for a radical simplification of 
the evolution equations before we can study the RG evolution of the Yukawa couplings in the MSSM with B and 
L violation. Motivated by the generational hierarchy of the conventional Higgs couplings, we shall assume that an 
analogous hierarchy amongst the different generations of B and L violating couplings exists. Thus, we shall retain 
only the couplings A233, A333, A233, and neglect the rest We note that B and L violating coupHngs to higher 
generations evolve more strongly because of larger Higgs couplings in their evolution equations, and hence could take 
larger values than the corresponding couplings to the lighter generations. We also note that the experimental upper 
limits are stronger for the B and L violating couplings with lower indices . With these assumptions wc can write 
the relevant renormalization group equations as P] 



167r^ 
167r2 
167r2 



d(-ln ^2) 

dhl 
d(-ln ^2) 

dhl 
d(-ln ^2) 



dX 



233 



d(-ln p?) 



16n^ rf^333 



d(~ln ^2) 

W)i"2 



d(-ln ^2) 



hUYs" + ^9i+^^9l~Gh't-hl-Xg 



333 2A2I3 ) I 



hi ( fgl + 35I + ^gl - hi - 6hl - hi 



6A333 2A2I3 I 5 



hi ( 3gl + -gl - Shi - Ahl - 4A233 - 3A'2 



333 ) I 



X2 

^^233 



3,9l + Tgl - K - 4A233 - 3A"2 



A333 ( y5l + 35: 



15 



5i 



2 hi 



233 



&hl - hi 



\II2 



•^233 6A333 2A233 I 1 



X"2 
^233 



8ff3 + 73i ~ '^h^ — — ^^333 ~ "^^233 



2A 



/2 



6A 



112 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
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where gi, 02,93 are the gauge couplings of U{1)y (in the GUT normaHzation), SU{2)l and SU{3)c gauge groups, 
respevtively, and /U is the running mass scale. We note that within the context of grand unified theories, one is led to 
the situation where baryon and lepton number violating Yukawa couplings may be related at the GUT scale, and one 
may no longer be able to set one or the other arbitrarily to zero. We, therefore, include both, the baryon and the lepton 
number violating couplings, in our RG equations. The evolution equations for the gauge couplings are not affected by 
the presence of B and L violating couplings at the one-loop level, and can be written, in the usual notation, as 

167^' :ir^^ = -^-9', * = 1, 2, 3, (11) 

d(— in 

with 

6, = (-, 1, -3). (12) 
The corresponding one-loop renormalization group equations for the gaugino masses M,;, z = 1, 2, 3 can be written as 



Defining 



IDTT^ lOTT^ 



\/2 \//2 

^ "16^' ^ ^16^' 

the solution of the RG equations (||) - ( p^ for the Yukawa and the B and L violating couplings can be written in a 
closed form [[121 



i+<>un-(o)/„*ft(c)<i(' 

where Mx is some large initial scale, and where Yk stands for the functions Yi {i — t, b, t), Y, Y' , and Y". Analogous 
notation holds for the functions Fk- The quantities Okk are the diagonal elements of the wave function anamolous 
dimension matrix, and are given by 

flfcfc = {6, 6, 4, 4, 6, 6}, (17) 
and the functions Fk are given by the set of integral equations 

Ft{t) = z r ^ J , (18) 

(1 + 6^(0) /„ Ft{t')dt'y/^l + 6Y'{0) /J F'{t')dt'y/^l + 6Y"{0) /J F"{t')dt'y/3 

Fbit) = . J . T ^ z z J , (19) 

(1 + 6Yt{0) /q Ft{t')dt'y/^1 + AY-riO) /g Fr{t')dt'y/*il + 6Y'{0) /J F'(t')di')'/'(1 + QY"iO) F"{t')dt'Yri 

FM) = : _ ^"'V z , (20) 

(1 + 6^(0) /q Fi,(t')dt'Y/''{l + 6F(0) /„ F{t')dt'){l + &Y'{0) /J F'{t')dVY/^ 

m = -~ 1 1 , (21) 

(1 + 4^(0) /J F^{t')dt'){l + 6F"(0) /J F"it')dt'y/2 



F'it) 



E'{t) 



(l + 6Yt(0) J^Ft{t')dt'Y/e(l + 6Yb{0) J^Ftit')dt')(l + 4Y.riO) F^{t')dt'Y/^{l + 49(0) F{t')dt'Y/* 
1 



(1 -f 6y"(o) p^F"{t')dt'Y^'-^ 



(22) 



F"{t) = . J . . J , (23) 

{l + GYt{Q) /jFt(i')dt')'/'(l + 6n(0) f^Ft{t')dt'Y/''{l + &Y'{Q) F'{t')dt'Y^--^ 
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where the functions Ek{t) (= Etit), Eb{t), Er{t), E{t),E'{t) and E"{t)) are given by 



with 



cu 



Ek{t)^Ylil + h&,{0)t) 



Cki/bi 



i=l 



5,(0) 



167r2 ■ 



i = l,2, 3, 



13 16\ / 7 16\ /9 

15' ^'Yj' ^"^^U'^'Tj' ^--U'^'"'' 



(24) 

(25) 
(26) 



3, 



7 16 

15"^' 3 



CA" 



0, 8 



(27) 



The solutions for the RG equations ( |Tll) for the gauge coupUngs and the gaugino masses (|T^ are well known and will 
not be repeated here. We note that (|16|) gives the exact solution for the Yukawa couplings, while Fk's in ( [l^ ) - ( |2^ ) 
should in principle be solved iteratively. 

In order to study the infrared fixed points for the Yukawa couplings, it is convenient to redefine 



Rt 



93 53 



233 



Rr — o : R — 9 

9i 53 



R'^ 



''333 

51 



^233 

9i 



(28) 



and retaining only the SU{3)c gauge coupling constant 173, we can write the renormalization group equations (^) 
( p^ ) for the Yukawa couplings as 



353 



35i 



Rt ~ ~ QRt — Rb — R' — 2i?" 
R}y — — Rt — 6i?f, — Rf — 6i?' — 2_R" 
Rr [-3 - 3i?6 - ARr -AR~ m'] , 



, dRt 

d9l 
, dRi, 

d9l 

^dRr 

3gl-^ = R [-3 - 4:Rr - 4i? - 3i?'] , 



Rt-GRb- Rr - R- 6i?' - 2i?" 



R" [5 - 2i?( - 2Rb - 2R' - 6i?"] 



(29) 
(30) 
(31) 
(32) 

(33) 
(34) 



There are no physically acceptable infrared fixed points for this set of RG equations with all the couplings attaining 
nontrivial values. The only nontrivial infrared fixed points are obtained by neglecting R' . Then, the true infrared 
fixed-points for these RG equations are M 



\ti , H^, Hi) 
{Rl, Rt) 
{R"\ Rl) 
{R"\ Rl) 



77 2 2 
102' 17' 17 
1 1 
3' 3 
19 1 
24' 8 
19 1 
24' 8 



r;^r* = (stable), 
, R*=R* = R"* = (unstable), 
, R*=R*=RI = (unstable), 
, R; = R* = R; ^0 (unstable). 



(35) 
(36) 
(37) 
(38) 



We note that the r Yukawa coupling Rr and the lepton number violating coupling R approach trivial fixed point values 
since there is no contribution from the SU{3) coupling 173 in their renormalization group equations. Also, because 
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of this reason the B- and L- violating couplings do not approach simultaneous non-trivial fixed points. Furthermore, 
the infrared fixed point values of the top- and bottom-quark Yukawa couphngs for the stable fixed point ( ^5|) are 
significantly different from the case when B and L is conserved Thus, the inclusion of B and L violation in 

MSSM has the effect of lowering the infrared fixed point values of the top- and bottom-quark Yukawa couplings. 
These are the important conclusions of our analysis. From (|3^) we note that the fixed point value for the top-quark 
Yukawa coupling translates into a top-quark (pole) mass of about mt ~ 70 sin /? GeV, which is incompatible with the 
measured value of the top mass, mt ~ 174 GeV, for any value of tan/3. Thus, the stable infrared fixed point ( |35| ) 
is not actually realised in nature. 

The infrared fixed points (IRFP's) that we have discussed above are the true IRFP's of the renormalization equa- 
tions for the Yukawa and baryon and lepton number violating couplings. However, these fixed points may not be 
reached in practice, the range between the large (GUT) scale and the weak scale being too small for the ratios to 
approach the fixed point values. In that case, the various Yukawa couplings may be determined by quasi-fixed point 
behaviour Q|, where the values of various couplings at the weak scale are independent of their values at the large scale, 
provided the Yukawa couplings at the large scale are large. More precisely, in the regime where the Yukawa couplings 
Yt{Q), Y(0), Yr(0), Y(0), y'jO), Y"(0) — > oo with their ratios fixed, it is legitimate to drop 1 in the denominators 
of the equations (Rq) and dlS) - (p3) so that the exact solutions for the Yukawa couplings approach the infrared 
quasi-fixed-point (IRQFP) defined by 



akkJoF^^''it')dt' 



(39) 



with 



Fr^'it) = 


Et{t) 


(40) 


(/o P^'^^i^'WY'^Uo F'Qpp{t')dt'y/(^{j* F"Qpp{t')dt'y/^ ' 


F^^it) = 


Ebit) 


(41) 


ilo Ft^''^it')dt'y/%J* F?^^{t')dt'Y/Hf^ F'QPP{t')dt'y/^{J* F"QPP{t')dt'Y/^ ' 


FQPP{t) = 


Er(t) 


(42) 


{f^F^^''{t')dt')y^UlFQPP{t')dt%jlF'QPP{t')dt'Y/^' 


FQPP{t) = 


Fit) 


(43) 


(/J F?^^{t')dt'){^l F"QPP{t')dt'Y/^ ' 


F'QPP^t) = 


E'it) 


(44) 


[f^F?^^{t')dt')y%f^F^^''{t')dt%f^F?^^{t')dt'^^^^^^ 


F"«^^(t) = 


E"{t) 


(45) 


UlFt{t')dt'Y'^UlFt{t')dt'Y/^U'^F'{t')dt'Y/^' 



We note that the result ( |39| ) follows immediately by dropping 1 from the denominator of (|lq), whereas the results 
( pO| ) - ( ^ follow from the corresponding equations ( p^ - ( |2^ ) by using an iterative procedure 15|. We stress 
here that both the dependence on the initial conditions for each Yukawa coupling as well as the dependence on the 
ratios of initial values of Yukawa couplings have completely dropped out of the runnings in Eqs.(|3^) and (^) - (|4^). 
In other words, the quasi-fixed-points ( p9|) are independent of whether the B and L violating couplings and the third 
generation Yukawa couplings are unified or not. The fact that the ratios of the various Yukawa couplings do not enter 
Eqs. (^) - (IH) implies that these results are valid for any tan/3 regime. 



B. Infrared fixed points for the trilinear soft supersymmetry breaking parameters 

We now consider the evolution equations for the soft supersymmetry breaking trilinear parameters in the potential 
(^. For these trilinear parameters we shall assume the same kind of generational hierarchy as was assumed for the 
corresponding Yukawa couplings. Thus, we shall consider only the highest generation trilinear coulings (^[7)33 = At, 
{Aohz = Ab, {Al)33 = Ar, (^a)233 = ^A, (^A')333 = Ay, (Aa")233 = Ay, Setting all Other elements equal to 
zero. As there is only one IRSFP ( p5| ) in the MSSM with B and L violation, we shall consider the IRFPs for the 
A parameters corresponding to this case only, i.e. for At, Ai,, A^, Ax and Ayi. Retaining only these parameters, 
and defining the ratios Ai = Ai/M^ {Ai — At, At, Ar, Ax, Ay), we can write the relevant renormalization group 
equations j|] for Ai as (neglecting the SU{2)l and U{1)y gauge couplings): 



6 



o 2d At 



35; 



,dAx 
' dgl 



16 



16 



(6i?t + 3)it - RbAb - 2R"A^ 



— - RtAt - i&Rb + 3)A - RrAr - 2R"A^ 



-3i?bib - [ARr + ^+ 3)i, - 



^i?,i,-(:^+4i?, + 3)iA 



2i?,tif - 2RbAi, - {6R" + 3)i> 



(46) 
(47) 
(48) 
(49) 
(50) 



One can obtain the exact solutions for RG equations ( p6| ) - ( pO[ ) for the trihnear parameters Ai analogous to the 
solutions (|6|) that we have obtained for the Yukawa couplings. The expressions for these solutions are lengthy, and 
will not be written down here. Instead we shall obtain approximate analytical solutions for these RG equations in 
Sec. Ill, and shall further study them numerically in Sec. IV. Nevertheless, we can obtain the true fixed points for the 
Ai parameters easily. These are 



{ai„,ai,a;,a;)^[i, 1, 1, 



(51) 



It is straightforward to show that this is the only infrared stable fixed point for the Ai paprameters. The stability 
of the fixed point ( |5T| ) also follows from the general results connecting the stability of a set of A parameters to the 
stability of the corresponding set of Yukawa couplings Q . We further note that the fixed point values for A), and At 
are the same as in MSSM with baryon number conservation [|l6|. However, the fixed point value for A^ is affected by 
the presence of the i?- violating parameter Ax" . 



III. ANALYTICAL SOLUTIONS 



Having obtained the infrared fixed points of the renormalization group equations for the Yukawa couplings and the 
soft supersymmetry breaking trilinear couplings, it is important to determine the approach to these fixed points. The 
rate of approach to the fixed points can be determined by solving the RG equations in the neighbourhood of the fixed 
points. Linearization of the RG equations (p9[)- ( p^ ) in the neighbourhood of the stable infrared fixed point ( |35| ) leads 
to the following approximate analytical solution for the Yukawa and B and L violating couplings: 



[0.5(i?to - Rbo) + 0.02i?,o] 



i?t - ^ - 0.002i?,o (4- 
17 V.930 

+ [0.48(i?to + Rbo) - 0.02i?^o - OMR'^ ~ 0.07] f 
+ [0.02(i?to + -Rbo) + 0.002i?^o + 0.06i?^' - 0.05] i^- 

\530 



at 

5io 



Rh 



17 



0.04i?^o 



at 

5io 



[0.5(i?to - Rw) + 0.02i?,o] 



at 

5io 



+ [0.48(i?to + i?6o) - 0.02i?^o - 0.06i?o - 0.07] 



(at 



-0.2 



+ [0.02(i?to + i?6o) + 0.002i?^o + 0.06i?o - 0.05] {^\ 

\530 / 

/ 2^-1 
V53O 



R = Ro 



al 
aio 



-0.2 



-1.18 



(52) 



-1.18 



(53) 
(54) 
(55) 
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R" = 0.75 - 0.005i?^o ^ 

.530 



2 N -1 

.91 



[-0.37(i?to + Rbo) + 0.02i?^o + 0.05< + 0.05] ( ^ 

.530 



2 \ -0-2 

9i 



o X -1.68 

55 



+ [0.37(i?to + i?6o) + 0.03i?^o + 0.95i?^' - 0.8] {^] , (56) 

V53o/ 

where Rto, RbOi RtO, Ro, R'o are the corresponding initial values at the large (UV) scale. These equations determine, 
in the leading order, how fast the infrared fixed point ( p5| ) is approached. The measure of infrared attraction is the 
size of the negative power in (g^/g^Q) in various terms in the expression for Ri. 

Similarly, liearization of the RG equations ( ^6| ) - ( ^0| ) for the trilinear couplings about the stable infrared fixed 
point ( ^l|) leads to the following approximate solution for these couplings: 

/ 2 \ -1-2 / 2 \ -2.57 

At = l + 0.5 (ito - Abo) ( ^ j +0.4 (iA"o - l) ( 

2 \ -1-21 

0.5 ' 



/ 2 \ — i.^i 

(l.2-(ito+Ao) + 0.8iA"o) (^^j , (57) 

/ 2 \ -1-2 / 2 \ -2.57 

ib = 1 - 0.5 (ito - ifco) j + 0.4 (Aa-o - l) ' 



530 



9 \ -1.21 

5i 



0.5(l.2-(Ato+Ao) + 0.8AA"o) ( ^ ) , (58) 

.530/ 



2 / s / „2 \ -1-2 , . / „2 N -2-57 



A. ^---0.3 (i.„ - i.„) I- + 0.03 (i.,„ - l) I- 



530/ ^ ^ \530 

0.47 + 0.03ito - 0.57Abo + ^ro + 0.19Aa"o 



2 \ -1 



do 



? \ -1.21 

53 



0.27 (1.2 - (Aio + A„) + 0.8Aa"o) ( ^ ) , (59) 

.530 . 



2 \ -1 
5,3 



Ax = Axo[^ , (60) 

.530 / 

9 \ -2.57 

55 



Ax" = 1 + ( /1a"o -1)1^ 

.530 



/ 2 \ -1.21 

+ 0.06(l.2-(i„ + Ao) + 0.8iA"o) (^^j , (61) 



where Ato,AbQ, Aro, Axq, Ayo are the initial values of the A parameters at the UV scale. 



IV. NUMERICAL RESULTS 



Having obtained analytically the true fixed points as well as the quasi-fixed points for the various couplings, as well 
as the approximate solutions to the corresponding RG equations, it is instructive to carry out a numerical study of 
these equations in order to illustrate the discussion of the previous sections. 

Since the RG equations are coupled differential equations with several fixed points, this implies that the system 
of equations (|2^) - (^) has an infrared fixed surface, the Rt — Rt — R" surface, Rt = R = 0. In principle, a 
three dimensional plot for the evolutions would reveal the infra-red fixed surfaces, with the sole infrared stable fixed 
point being approached by the fiow from different directions. However, such a demonstration is rather cumbersome. 
We shall instead concentrate on the two dimensional projections of this three dimensional surface to illustrate the 
approach to the fixed point. 

For the purposes of our numerical analysis, we set the unification scale Mx — 10^^ GeV, the unified gauge coupling 
ctciMx) = 1/24.5 and a effective supersymmetry scale of Ms = 1 TeV. We shall use the one- loop RG equations in 
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the numerical analysis, since the difference between the one-loop and two-loop results for the infrared fixed points |l^] 
is less than 10%. Reinserting the electroweak gauge couplings 51,52 in the RG equations ( p9|) ~ (^, the fixed points 
are no longer found to be exact. In order to determine the approximate positions of the fixed points, we shall maintain 
running g2(/i) and follow the prescription that the approximate infrared fixed point in the Rt — Rb — R" space 

is determined as the unique point in the evolution path from /i — Mx to fi — Ms which has the same value at 
fi = Mx and at /i = Ms- The results depend somewhat on the evolution path, and in particular on the scale Mx, the 
deviations from the values in (|3^) being larger, the larger the scale Mx- We note that the idea of grand unification 
is not important for our results. The infrared fixed point and the infrared fixed surface are the properties of minimal 
supersymmetric standard model with baryon and lepton number violation. Note that it is only in that they are not 
exact but approximate do they depend on Mx- 

Our main objective here is to study the renormalization group flow in the coupling ratios Rt, Rb, R" from Mx to 
Ms, to display the fixed surfaces (projections of three dimensional surface onto two dimensions), and to locate the 
infrared fixed point. We start with ht and hb each taking the values (9, 1, 0.5) to ensure a large range in the variation 
of these parameters and A233 fixed at a reasonably large value of 1.1 for the computations that we perform. We 
evolve the Yukawa couplings down to the IR scale using the one-loop RG equations and evaluate the coupling ratios 
Rt, Rb, R" - The values of the couplings so obtained at the IR scale, viz., Ms, are then used as the input at the UV 
scale, Mx, and evolved down to the IR scale and then iterated again. This procedure is repeated until all the final 
points meet in the infrared fixed point. The results of these iterations are then projected onto the Rt — Rb, Rt — R" 
and Rb — R" planes, respectively. The results so obtained are presented in Figs. 1(a), 1(b) and 1(c), respectively. It 
may be observed that the approach to the fixed point of the ratio R" is rather rapid since at the outset this quantity 
was chosen close to it. The approach of the ratios Rt and Rb is significantly more dramatic since they are required to 
approach the fixed point from rather disparate values. This clearly illustrates the strong attraction to the sole IRSFP. 

The end result of the series of iterations that we have performed lead to the following infrared fixed point values 
for the coupling ratios: 

R"Q^P ~ 0.65, 
i?^^^ ~ 0.59, 

Rf^^ ~ 0.59. (62) 

This infrared fixed point and the fixed surfaces are as expected from the preceeding analytical considerations with, 
however, shifted position of the fixed point due to the effect of including the electroweak gauge couplings in our 
numerical iterations. It is, in fact, the quasi fixed point for the respective coupling ratios. The fixed point ( |62|) implies 
the approximate top-bottom Yukawa unification at all scales /i. 

The fixed point values in ( |62[ ) translate into the fixed point values for the top-, bottom-quark Yukawa couplings, 
and the baryon number violating coupling: 

^sr - 0.90, 
h^^^ ~ 0.86, 

hf^^ ~ 0.86. (63) 

These quasi-fixed-point values for the Yukawa couplings are not significantly different from those obtained in a 
situation when the r Yukawa coupling was ignored, and when the baryon and lepton number violating couplings were 
considered separately in the fixed point analysis . Since the quasi-fixed points are reached for large initial values of 
the couplings at the GUT scale, these reflect on the assumption of purturbative unitarity, or the absence of Landau 
poles, of the corresponding couplings. The quasi- flxed points (|63|), therefore, provide an upper bound on the relevant 

Yukawa and the baryon and lepton number violating couplings. From our analysis we, thus, conclude that A233 ~ 1 
in a process independent manner. 

We have carried an analogous numerical study for the soft supersymmetry breaking trilinear couplings. Here the 
choice of inputs is larger and we economize on the possibilities, however, ensuring that the main focus of the analysis 
is not lost. In Figs. 2(a), 2(b) and 2(c) we illustrate the infrared attraction for the trilinear couplings with a choice of 
Yukawa and the baryon number violating couplings hb = A233 = 1-1 and ht = 1.65 (in order to break the near isospin 
invariance of the top- and bttom-quark couplings), and universal boundary conditions for the trilinear couplings 
±(3,2, l,0.5)Afi/2, 0, with M1/2 the universal gaugino mass, for the eight different computations that we perform. 
The results for the Yukawa couplings and the trilinear couplings at the infrared scale are then fed back to the RG 
equations at the UV scale and then evolved again. We exhibit the resulting attraction in terms of the projections 
onto the At — Ab, At — Ayi and Ab — A\n planes in Fig. 2. The dramatic focussing property of the renormalization 
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group equations for the A parameters is clearly observed in these results. The results of this iteration process results 
in the following predictions for the quasi-fixed points for the A parameters: 

A'^,^^ ~ 0.95 mg, 
Af^^ ~ 0.66 mg, 

A^^^~0.67mg, (64) 
where rrig is the gluino mass (= M3) at the weak scale. These quasi-fixed point values for the A parameters must 



be compared with the true fixed point values (51). We note that the q uasi-fixed-point values (^4|) provide a lower 
bound on the corresponding A parameters, whereas the true fixed-point values (|5l| ) represent an upper bound on 
these parameters. From this analysis we are able to constrain these A parameters in the following model independent 
manner: 

rUg 

m~g 

0.67 - — ^1, (65) 



V. SUMMARY AND CONCLUSIONS 



We have carried out a detailed study of the infrared fixed point structure of the minimal supersymmetric standard 
model with the third generation Yukawa couplings and with highest generation baryon and lepton number violation. 
We have obtained the infrared fixed points for such a model, and shown that there is no physically acceptable infrared 
fixed point with both baryon and lepton number violating couplings approaching a nontrivial fixed point. The 
simultaneous nontrivial fixed point for the top- and bottom-quark Yukawa couplings, and the i?-violating coupling 
A233 is the only true fixed point that is stable in the infrared region. We have obtained approximate analytical 
solutions to the RG equations near the stable fixed point, which illustrate the approach to the fixed point. We have 
also derived the exact solutions of the RG equations of such a model in a closed form, from which we have obtained 
infrared quasi-fixed point solutions for the various couplings. The quasi-fixed-points are realised at the weak scale 
when the initial couplings at the GUT scale are large. These fixed points, thus, reflect on the assumption of the 
perturbative unitarity of the corresponding couplings. We have carried out the corresponding RG analysis of the 
trilinear soft supersymmetry breaking couplings, and obtained the infrared stable fixed point for these couplings, and 
the approach to the fixed point analytically. 

Since the true fixed points may not be reached at the electroweak scale, we have also studied the numerical 
solutions of the RG equations, and obtained the infrared fixed surfaces, and demonstrated the convergence of the RG 
fiow towards the fixed point. From this analysis we have obtained a process independent upper bound on the highest 
generation baryon number violating coupling, A233 ^ 1. 

Our study of the RG flow of the Yukawa couplings has been complemented by the corresponding numerical study 
of the RG equations for the soft supersymmetry breaking trilinear couplings, and the demonstration of the rapid 
convergence towards the fixed point for these couplings. From this analysis, we have constrained the A parameters to 
be A\ii/mg ~ 1, 0.66 ~ At/rrig ^ 1, 0.67 ~ Ab/rrig ^ 1. We emphasize that our results are independent of whether 
or not there is a unification of the couplings at some large scale or not. The infrared fixed points and the infrared 
fixed surfaces are the property of the MSSM with baryon and lepton number violation. Since they are approximate, 
only their actual value depends on the large scale. 
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Figure Captions 

Fig. 1 (a) The values of the ratios Rb and Rt at the end of first four iterations for the inputs ht = 9, 1,0.5, hb = 

9,1,0.5, and X'^^^ = 1.1. 

Fig. 1 (b) The values of the ratios R" and Rt at the end of first four iterations for the inputs ht = 9, 1,0.5, hb = 
9,1,0.5, and A^'33 = 1.1. 

Fig. 1 (c) The values of the ratios R" and Rb at the end of first four iterations for the inputs ht = 9, 1,0.5, hb = 
9,1,0.5, and A^'33 = 1.1. 

Fig. 2 (a) The values of the ratios Ab and Ai at the end of first four iterations for the inputs hb = A233 = 1.1 and 
ht = 1.65, and universal trilinear couplings A — ±(3, 2, 1, 0.5)Mi/2, 0. 

Fig. 2 (b) The values of the ratios A\» and At at the end of first four iterations for the inputs hb = A233 = 1.1 and 
ht = 1.65, and universal trihnear couplings A = ±(3, 2, 1, 0.5)Mi/27 0. 

Fig. 2 (c) The values of the ratios Ay and Ab at the end of first four iterations for the inputs hb = A233 = 1.1 and 
ht = 1.65, and universal trilinear couplings A = ±(3,2, l,0.5)Mi/2, 0. 
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